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ABSTRACT 


Defining  the  averaging  time  required  for  measuring  meaningful  turbulence  statistics  is  a  central  problem 
in  boundary-layer  meteorology.  Path-averaging  scintillation  instruments  are  presumed  to  confer  some 
time-averaging  benefits  when  the  objective  is  to  measure  surface  fluxes,  but  that  hypothesis  has  not  been 
tested  definitively.  This  study  uses  scintillometer  measurements  of  the  inner  scale  of  turbulence  {(■  ()  and 
the  reflective  index  structure  parameter  (C^)  collected  during  SHEBA  (the  experiment  to  study  the  Sur¬ 
face  Heat  Budget  of  the  Arctic  Ocean)  to  investigate  this  question  of  required  averaging  time.  The  first 
conclusion  is  that  the  beta  probability  distribution  is  useful  for  representing  €„  and  measurements. 
Consequently,  beta  distributions  are  used  to  set  confidence  limits  on  C  „  and  I  o  values  obtained  over 
various  averaging  periods.  When  flie  Cn  and  I  q  time  series  are  stationary,  a  short-term  average  of  Cn  or 
1 0  can  be  as  accurate  as  a  long-term  average.  But,  as  with  point  measurements,  when  time  series  of  path- 
averaged  Cn  or^ovalues  are  nonstationary,  turbulent  surface  fluxes  inferred  firom  these  Cn  and  .^o^alues 
can  be  variable  and  uncertain — problems  that  path-averaging  was  presumed  to  mitigate.  Since 
nonstationarity  turns  out  to  be  a  limiting  condition,  the  last  topic  is  quantifying  tiie  nonstationarity  with  a 
published  nonstationarity  ratio  and  also  by  simply  counting  zero-crossings  in  the  time  series. 
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P.  OLA  G.  PERSSON,  AND  PETER  S.  GUEST 


1  INTRODUCTION 

The  wide-ranging  fields  of  atmospheric  optics  and  boimdary-layer  meteor¬ 
ology  intersect  when  the  subject  is  optical  turbulence.  The  classic  works  by 
Tatarskii  (1961, 1971)  and  Chernov  (1967)  define  the  scope  of  this  intersection. 
More  recently,  Andreas  (1990)  assembled  classic  papers  from  both  disciplines 
that  focus  on  obtaining  turbulent  fluxes  in  the  atmospheric  boundary  l^er  from 
the  propagation  statistics  obtained  at  optical  or  other  electromagnetic  (EM) 
wavelengths. 

This  idea  of  using  the  propagation  statistics  of  EM  waves — especially  those 
that  quantify  scintillation — ^to  infer  turbulent  properties  of  the  atmospheric 
boundary  l^er  has  been  around  for  about  35  years  (e.g.,  Strohbehn  1970,  Gray 
and  Waterman  1970,  Gurvich  et  al.  1974).  Because  of  the  need  to  evaluate  the 
surface  energy  budget  over  terrain  that  m^  be  inhomogeneous,  the  emphasis 
in  the  last  20  years  has  been  on  using  scintillation  statistics  to  infer  the  surface 
fluxes  of  heat,  moisture,  and  momentum  (e.g.,  Kohsiek  and  Herben  1983,  Hill  et 
al.  1988,  1992a,  1992b,  Andreas  1989b,  1992,  Thiermarm  and  Grassl  1992, 
Green  et  al.  1994,  De  Bruin  et  al.  1995,  Nieveen  and  Green  1999). 

Two  propagation  statistics  are  commonly  used  to  estimate  turbulent  fluxes  in 
the  atmospheric  surface  l^er:  the  refractive  index  structure  parameter,  €„ ,  and 
the  inner  scale  of  turbulence,  •  Th®  refractive  index  structure  parameter  sets 
the  level  of  the  inertial-convective  subrange  for  the  one-dimensional  refractive 
index  spectrum,  On  (e.g.,  Andreas  1987): 


<D„(k)  =  0.249CJk-'''  , 


(1.1) 


2 


ERDC/CRRELTR-03-24 


where  k  is  the  magnitude  of  the  turbulence  wavemunber.  The  refractive  index 
spectrum  in  the  inertial-convective  subrange,  however,  can  also  be  expressed  in 
more  familiar  meteorological  variables: 

0„(k)  =  P„N„s-‘'^k-’'\  (1.2) 

where  Nn  is  the  dissipation  rate  of  the  variance  in  refractive  index,  e  is  the  dissi¬ 
pation  rate  of  turbulence  kinetic  energy,  and  Pn  (~  0.4,  Andreas  1987)  is  the 
Kolmogorov  (or  Corrsin)  constant  for  the  refractive  index  spectrum. 

For  the  optical  wavelengths  that  we  use,  relates  directly  to  the  tempera¬ 
ture  structure  parameter  Cj  through 

=  A^(?.,P,T,Q)(4  .  (1.3) 

Here,  A  is  a  known  function  of  the  EM  wavelength  (1)  and  the  average  baro¬ 
metric  pressure  (P),  temperature  (T),  and  specific  humidity  (Q)  (e.g.,  Andreas 
1988b). 

In  (1.3),  we  ignore  contributions  from  the  humidity  and  temperature- 
humidity  structure  parameters.  These  contributions  are  usually  negligible  unless 
the  absolute  value  of  the  Bowen  ratio  is  small  (e.g.,  Wesely  and  Alcarez  1973, 
Wesely  1976,  Thiermann  and  Grassl  1992,  Hill  et  al.  1992b).  Over  sea  ice,  the 
magnitude  of  the  Bowen  ratio  is  almost  always  large  enough  for  (1 .3)  to  be  a 
good  approximation  (Andreas  and  Cash  1996). 

Because  the  temperature  spectrum  also  obeys  equations  like  (1.1)  and  (1.2), 
(1.3)  can  be  written  as 

Cj  =  A'(Pt/0.249)Nt8-‘'’,  (1.4) 

where  Pt  is  the  Kolmogorov  constant  for  the  temperature  spectrum,  and  Nx  is  the 
dissipation  rate  of  temperature  variance.  As  a  result,  is  related  to  dissipation 
rates  that  are  commonly  used  in  the  so-called  inertial-dissipation  method  for  flux 
estimation  (e.g.,  Fairall  and  Larsen  1986). 

The  inner  scale  of  turbulence  likewise  relates  to  features  of  the  refractive 
index  spectrum;  £q  is  approximately  the  turbulent  eddy  size  that  separates  the 
inertial-convective  subrange  from  die  dissipation  region  and  is  roughly  seven 
times  larger  than  the  Kolmogorov  microscale.  Because  at  optical  wavelengths 
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the  temperature  spectrum  and  the  refractive  index  spectrum  have  identical  shapes 
at  large  wavenumbers,  we  can  write  (Hill  and  Clifford  1978) 


^0  =  [9r(l/3)pTDf' (1.5) 

where  F  is  the  gamma  function  and  D  is  the  thermal  diffusivity  of  air.  Like 
Cn ,  the  propagation  statistic  is  related  to  a  dissipation  value  (i.e.,  s)  that  is 

required  for  using  the  inertial-dissipation  method  to  estimate  the  turbulent  surface 
stress. 

Two  routes  are  available  now  for  estimating  the  surface  stress  (or  momentum 
flux)  and  the  sensible  heat  flux  from  fliese  scintillation  measurements.  One  way 
is  to  use  the  dissipation  rates  e  and  Nt  computed  from  (1.5)  and  (1.4),  respective¬ 
ly,  in  the  traditional  inertial-dissipation  equations  (e.g.,  Large  and  Pond  1982, 
Fairall  and  Larsen  1986,  Edson  et  al.  1991).  The  second  w^  is  to  use  the 
measmed  refractive  index  structure  parameter  in  a  similarity  relation  for  the 
temperature  structure  parameter. 

For  this  latter  method,  the  two  relevant  equations  are  (e.g.,  Andreas  1988b, 
Thiermaim  and  Grassl  1992,  Hill  1997) 


and 


t.  - 


,2/3, 


-il/2 


A^g(C) 


u*  = 


&kz 


,1/3 


(1.6) 


(1.7) 


Here,  z  is  the  path  height,  k  is  the  von  Karman  constant,  and  g(C)  and  (|)e(C)  are 
semi-empirical  Monin-Obukhov  similarity  functions  of  the  stability  parameter 
^  =  z/L,  where  L  is  the  Obukhov  length  (see  Appendix  A).  Also  in  these,  u,  is 
the  friction  velocity  and  t*  is  a  temperature  flux  scale  such  that  T  =  p  u,  is  the 
surface  stress  and  =  -  pCp  u,  t,  is  the  smface  sensible  heat  flux,  where  p  is 

the  air  density  and  Cp  is  the  specific  heat  of  air  at  constant  pressure.  Note  fliat, 
because  of  the  square  root,  (1.6)  is  ambiguous  as  to  the  sign  of  t. ;  we  must 
evaluate  this  by  other  means. 
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Using  scintillation  statistics  to  estimate  the  turbulent  surface  fluxes  is  attrac¬ 
tive  because  the  EM  waves  propagate  over  some  finite  path.  The  propagation 
statistics  are  thus  path-averaged,  as  would  be  any  turbulence  quantity  derived 
from  them.  Presumably,  sueh  path-averaging  would  mitigate  the  effects  of  minor 
nonstationarity  or  surface  heterogeneity  that  often  confound  turbulence  meas¬ 
urements  made  with  point  sensors  (Schmid  and  Oke  1990,  Horst  and  Weil  1992, 
1994). 

A  seeond  supposed  benefit  of  using  EM  propagation  measurements  of 
and  £q  to  determine  u»  and  t»  is  that  statistieally  reliable  values  of  and  £q 
can  be  obtained  with  veiy  short  averaging  times,  on  die  order  of  one  minute  or 
less  (Wyngaard  and  Clifford  1978,  Frehlich  1988).  By  extension,  this  result 
would  seem  to  suggest  that  x  and  Hs  could  also  be  measured  using  path-averaging 
instruments  with  much  shorter  averaging  times  than  the  30-60  minutes  of  aver¬ 
aging  typieally  required  for  point  measurements  (Haugen  et  al.  1971,  Wyngaard 
1973,  Sreenivasan  et  al.  1978). 

The  rationale  for  this  idea  that  short  flux-averaging  times  are  possible  is  that 
a  path-averaging  sensor  samples  many  more  turbulent  eddies  per  unit  time  than  a 
point  sensor  (Andreas  1988a).  Some  of  those  eddies  would,  presumably,  be  the 
larger  ones  that  manifest  as  nonstationarity  or  intermittency  in  a  time  series  of 
point  measurements.  The  path-averaging  supposedly  would  smooth  out  the 
effects  of  these  large  eddies  and  thereby  reduce  the  sampling  error  common  in 
point  measurements  of  turbulent  fluxes. 

Even  if  this  scenario  is  accurate,  though,  we  are  still  skeptical  that  path- 
averaged  scintillation  statistics  can  routinely  yield  meaningful  estimates  of  x  and 
Hs  for  averaging  times  on  the  order  of  minutes.  Obtaining  these  turbulent  flux 
estimates  requires  convolving  the  scintillation  statistics  and  £q  (or,  equiva¬ 
lently,  e)  iteratively  with  the  Monin-Obukhov  similarity  fimctions  g(C)  3nd 
(|)g  (i^)  in  (1.6)  and  (1.7).  All  evaluations  of  these  similarity  fimetions  that  we 
know  of,  however,  are  based  on  point  measurements  of  mean  meteorological  and 
turbulence  quantities  averaged  for  30-60  minutes  (e.g.,  Businger  et  al.  1971, 

Dyer  and  Bradley  1982,  Hogstrom  1988,  Oncley  et  al.  1996,  Edson  and  Fairall 
1998).  It  is  not  obvious  to  us  that,  in  general,  combining  minute  averages  of 
and  £(,  with  functions  based  on  at  least  30-minute  averages  yields  meaningful 
values  of  x  and  Hs  (cf.  Andreas  1988a).  If  short-term  averages  of  €„  and  £q 
are  representative  of  time  series  averaged  for,  say,  30-60  minutes,  however, 
using  these  short-term  averages  in  (1.6)  and  (1.7)  to  estimate  x  and  Hs  may  be 
defensible. 

Here  we  address  these  issues  of  how  to  properly  average  scintillation-derived 
values  of  and  £q  for  the  purpose  of  estimating  the  turbulent  surface  fluxes. 
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First  we  review  the  time-averaging  constraints  for  evaluating  Monin-Obukhov 
similarity  functions  to  highlight  our  concern  over  combining  similarity  fimctions 
with  short-term  scintillometer  averages.  Next  we  show  seasonal  histograms  for 
and  values  based  on  600  hours  of  scintillometer  data  collected  during 

SHEBA,  the  experiment  to  study  the  Surface  Heat  Budget  of  the  Arctic  Ocean. 
We  find  that  beta  probability  distributions  fit  both  variables  well.  We  also  show 
that  beta  distributions  fit  hour-long  series  of  and  values  that  contain  only 
60  points.  We  therefore  use  beta  distributions  to  compute  error  bars  for  several 
representative  hour-long  time  series  of  SHEBA  and  values.  This  error 
analysis  reveals  that  intermittency  and  nonstationarity  can  have  pronounced 
effects  on  these  path-averaged  €„  and  £q  values,  contrary  to  some  of  our  earlier 
speculation.  As  a  result,  we  reiterate  that  short-term  averages  of  €„  and  £q  can 
yield  reliable  short-term  averages  of  x  and  Hs  only  for  quasi-stationary  time 


senes. 
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2  OBSERVATIONS 

SHEBA  was  a  multidisciplinary,  yearlong  experiment  on  drifting  sea  ice  in 
the  Beaufort  Gyre  (Perovich  et  al.  1999,  Uttal  et  al.  2002).  Andreas  et  al.  (1999) 
and  Persson  et  al.  (2002)  describe  our  SHEBA  program.  Andreas  et  al.  (2000) 
make  a  preliminary  report  of  our  current  analysis,  and  Andreas  et  al.  (2003)  give 
a  shorter  version  of  this  report. 

Our  SHEBA  scintillometer  (Fig.  1)  was  an  SLS20  system  made  by  Scintec 
Atmospharemnesstechnik  GmbH  of  Tubingen,  Germany  (Thiermaim  1992).  Its 
source  is  a  laser  of  0.685-pm  wavelength.  We  operated  this  system  over  a  350-m 
path  at  a  height  of  2.88  m  for  our  1997  SHEBA  measurements  and  over  a  300-m 
path  at  a  height  of  2.60  m  for  our  1998  measurements.  The  surface  was  generally 
snow-covered  sea  ice. 


Figure  1.  The  scintillometer  receiver  at  the  SHEBA  ice  camp.  The  scintil¬ 
lometer  transmitter  is  in  the  distance. 


At  the  scintillometer  source,  the  laser  beam  is  split  into  two  beams  with  their 
centers  separated  by  2.7  mm.  The  scintillometer  receiver  thus  has  two  detectors, 
also  separated  by  2.7  mm.  The  correlation  between  the  intensity  fluctuations  at 


Flux  Averaging 


7 


the  two  detectors  is  a  measure  of  .  The  intensity  fluctuations  measimed  by 
either  detector  are  related  to  . 

The  fundamental  data  that  the  SLS20  reports  are  minute  averages  of  €„  and 
and  a  data-quality  number  designated  NOK  (for  number  okay).  To  obtain 
these  averages  and  NOK,  the  system  software  divides  each  minute  into  ten  6-s 
blocks.  During  the  first  block  of  each  minute,  the  software  turns  the  laser  off,  and 
the  receiver  measures  the  background.  The  system  software  then  applies  this 
background  calibration  to  the  nine  subsequent  6-s  blocks  in  that  minute.  Each  of 
these  nine  blocks  yields  an  individual  measurement  of  €„  and  ;  the  software 
also  evaluates  the  quality  of  each  6-s  block.  If  a  block  passes  the  quality-control 
criteria,  its  €„  and  £q  values  are  used  for  computing  the  minute  average;  other¬ 
wise,  that  block’s  values  are  ignored.  NOK  reports  how  many  of  the  nine  avail¬ 
able  data  blocks  during  the  minute  were  used  to  compute  the  minute  averages  of 
and  .f  0 .  To  compile  the  histograms  that  we  show  in  Section  4,  we  rejected 
any  minute-averaged  €„  and  values  that  were  not  based  on  at  least  three 
valid  6-s  blocks  during  that  minute.  That  is,  we  retained  the  minute-averaged 
and  £q  values  for  further  analysis  only  if  NOK,  which  is  reported  as  a  percent¬ 
age,  was  at  least  33%  [=  (3/9)xl00%]. 

Figure  2  shows  a  daylong  time  series  of  typical  minute-averaged  €„  and  £q 
values.  Generally,  £q  is  inversely  related  to  the  wind  speed;  higher  wind  speeds 
increase  the  dissipation  rate  of  turbulent  kinetic  energy,  which,  in  turn,  is  asso¬ 
ciated  with  smaller  £q  values  [see  (1.5)].  The  behavior  of  is  not  as  easy  to 

understand.  On  the  one  hand,  higher  wind  speeds  usually  foster  mixing,  which 
erodes  any  vertical  temperature  gradient,  so  decreases  (cf.  Frederickson  et  al. 
2000).  But  on  the  other  hand,  if  some  process  such  as  radiative  heating  or  cooling 
can  maintain  a  vertical  temperature  gradient,  a  higher  wind  speed  would  increase 
the  magnitude  of  the  sensible  heat  flux  and,  thereby,  increase  C„ . 

Later  we  will  compute  statistics  from  series  of  and  £q  such  as  in  Figure 
2  and  derive  confidence  limits  from  these  statistics.  Such  analyses  get  compli¬ 
cated  if  adjacent  samples  in  the  time  series  are  correlated.  The  autocorrelation 
function  is  typically  used  to  quantify  the  decorrelation  time  of  a  random  process. 
The  integral  of  the  autocorrelation  function  over  all  time  lags  defines  a  time  3 
called  the  integral  scale  (e.g.,  Lumley  and  Panofsky  1964,  p.  36£)  that  quantifies 
this  decorrelation  time.  Sample  values  separated  in  time  by  intervals  of  3  or 
longer  are  essentially  uncorrelated  or  independent. 

As  an  estimate  for  the  integral  scale  in  the  neutrally  stratified  atmospheric 
siuface  layer,  Wyngaard  (1973)  used 
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Figure  2.  A  day  of  minute-averaged  inner  scale  {t^)  and  refractive  index 
structure  parameter  ( )  values  at  a  height  of  3  m  from  SHEBA. 

From  high-frequency  time  series  of  temperature,  hmnidity,  and  velocity 
components,  Sreenivasan  et  al.  (1978)  actually  computed  autocorrelation  and 
cross-correlation  functions.  From  the  autocorrelation  functions  for  the  individual 
variables,  they  deduced 
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3  =  —.  (2.3) 

U 

Their  analysis  of  the  cross-correlations,  on  the  other  hand,  yielded  a  multipli¬ 
cative  constant  in  (2.3)  that  was  typically  1.2. 

To  be  conservative,  we  will  use  (2.3)  to  estimate  the  integral  scale.  Before 
considering  our  scintillometer  data,  though,  we  can  confirm  the  usefulness  of 
(2.3)  with  computations  reported  by  Trevino  and  Andreas  (2000).  For  three  short 
records  of  the  longitudinal  velocity  component  from  a  sonic  anemometer 
sampling  at  10  Hz,  they  evaluated  a  quantity  they  called  the  “memory,”  which 
approximates  the  decorrelation  time  or  3.  These  memory  values  ranged  from 
0.88  to  1.46  s.  For  their  measmements,  z  was  4  m  and  U  was  about  2  m  s“\ 
Equation  (2.3)  thus  predicts  3  =  8  s  for  these  data,  which  is  a  conservative 
estimate  of  the  computed  decorrelation  times. 

Our  scintillometer  operated  at  heights  a  bit  under  3  m;  in  (2.3),  we  will  thus 
use  3  m  for  z.  The  integral  scale  clearly  increases  as  the  wind  speed  decreases; 
hence,  the  largest  estimate  of  the  decorrelation  time  is  associated  with  the  lightest 
winds.  In  fewer  than  10  of  the  600  hours  of  scintillometer  data  that  survived  our 
quality  controls  was  the  average  hourly  wind  speed  less  than  1ms"’.  Substituting 
1  m  s"'  and  z  =  3  m  in  (2.3)  gives  12  s  as  a  conservative  estimate  for  the  decor¬ 
relation  time  of  our  raw  scintillometer  data.  For  most  of  the  data,  the  decorrela¬ 
tion  time  is  much  shorter  because  the  wind  speed  was  higher.  Since  the  only 
scintillometer  data  we  use  in  our  subsequent  analyses  are  minute  averages,  we 
can  be  confident  that  these  minute-averaged  and  £q  values  are  independent. 
We  thus  need  not  worry  about  the  complicating  effects  of  correlated  data  in  our 
statistical  analysis. 
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3  AVERAGING  TO  ESTABLISH  MONIN-OBUKHOV 
SIMILARITY  THEORY 

Equations  (1.6)  and  (1.7),  which  are  needed  to  convert  scintillometer 
measurements  to  flux  estimates,  derive  from  Monin-Obukhov  similarity  theory. 
Evaluating  the  similarity  functions  g(Q  and  <|)e(Q  that  make  (1.6)  and  (1.7)  work 
requires  measuring  both  means  and  covariances.  In  particular,  we  must  measure 
the  mean  wind  speed  (U)  and  temperature  (T)  profiles  and  the  eovariances  uw 
( = -  uj )  and  wt  ( = -  u»  t*),  where  die  overbar  denotes  averaging.  To  have 
confidence  in  the  resulting  similarity  fimctions,  we  must  have  eonfidence  in  the 
statistical  properties  of  these  averages.  That  is,  we  must  average  long  enough  to 
ensure  reliable  statistics. 

Lumley  and  Panofsky  (1964,  p.  35ff.),  Wyngaard  (1973),  Sreenivasan  et  al. 
(1978),  Andreas  (1988a),  and  Lensehow  et  al.  (1994),  among  others,  estimate 
how  low  is  long  enough  to  average.  Here  we  rely  on  the  analyses  by  Wyngaard 
and  Sreenivasan  et  al. 

Building  on  work  by  Lumley  and  Panofsky  (1964),  Wyngaard  (1973) 
^timated^the  averaging  time  Tx  required  to  reduce  the  error  in  a  measured  mean 
X  to  X .  He  obtained 


(3.1) 


where  3x  denotes  the  integral  scale  of  variable  x,  and  x*  is  the  ensemble 
variance  about  x .  As  we  mentioned  earlier,  Sreenivasan  et  al.  (1978)  evaluated 
5  for  quantities  like  U,  T,  uw ,  and  wt .  In  general,  they  foxmd 


(3.2) 


where  otx  is  3.9, 4.9, 1.2,  and  1.2  for  U,  T,  uw ,  and  wt ,  respectively. 

Evaluating  the  Monin-Obukhov  similarity  fimctions  requires  measuring 
the  vertical  gradients  in  wind  speed  and  potential  temperature.  We  must,  thus, 
require  an  accuracy  in  individual  wind  speed  and  temperature  measmements 
of,  say,  5  cm  s”’  and  0.1  K,  respectively.  For  a  5-m  s"*  wind  speed,  this  makes 
6u  =  0.01;  for  an  average  temperature  of  293  K,  6t  =  0.00034. 
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To  estimate  the  averaging  time  for  finding  the  vertical  gradient  in  wind 
speed,  combine  (3.1)  and  (3.2)  to  get 


Tu 


=  2x10'' 


p.9z^ 

(ol] 

u  J 

tu’J 

(3.3) 


where  Oy  is  the  variance  in  the  longitudinal  component  of  the  wind  vector. 
Typically,  Ou/U  ~  0.1  (e.g.,  Arya  1988).  Hence,  for  a  10-m  profiling  tower 
in  a  5-m  s”'  wind,  Tu  «  1600  s  «  26  minutes  is  an  approximate  averaging  time. 

Similarly,  to  measure  the  mean  vertical  temperature  gradient,  we  must 
average  according  to 


Tj  =  1.7x10’ 


f  4.9z') 

[  u  J 

T2 

J 

(3.4) 


When  the  sensible  heat  flux  is  on  the  order  of  100  W  m"^,  is  about  0.5  K^. 
Hence,  for  an  average  temperature  of  293  K  and  the  other  conditions  as  above, 
Tt  ~  970  s  ~  16  minutes  is  an  approximate  averaging  time. 

That  is,  measuring  the  mean  profiles  necessary  to  determine  the  Monin- 
Obukhov  similarity  functions  requires  averaging  for  15  to  30  minutes. 

To  estimate  the  flux-averaging  time,  Wyngaard  (1973)  and  Sreenivasan  et  al. 
(1978)  wrote  (3.1)  as 


Cu 


X  v2  2 

(wxj  -  wx 

=2  ’ 
WX 


(3.5) 


where  x  is  now  the  turbulent  fluctuation  in  eifiier  the  longitudinal  velocity  (u)  or 
temperature  (t).  From  Kansas  data,  Wyngaard  estimated  the  bracketed  quantity  in 
(3.5)  to  be  about  10  for  both  uw  and  wt  for  near-neutral  stratification.  Sreeni¬ 
vasan  et  al.  report  that  this  quantity  is  15  and  32,  respectively,  for  uw  and  wt . 

An  error  of  10%  (i.e.,  6wx  =  0.1)  is  about  the  best  that  has  been  demonstrated 
for  eddy-correlation  measurements  of  the  turbulent  surface  fluxes.  Thus,  for  a 
measurement  height  of  5  m,  a  mean  Avind  speed  qfS  m  s"',  and  with  =1.2 
from  Sreenivasan  et  al.  (1978)  for  both  uw  and  wt ,  (3.5)  gives  T  UW  Tvrt  ®  40 
minutes  when  we  use  Wyngaard’s  (1973)  value  for  the  bracketed  term  in  (3.5). 
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When  we  use  the  estimates  for  this  term  from  Sreenivasan  et  al., 

Tuw  »  60  minutes  and  T^t »  128  minutes. 

In  summary,  using  eddy  correlation  to  measure  the  momentum  and  heat 
fluxes  to  within  10%  for  evaluating  the  Monin-Obukhov  similarity  functions 
requires  roughly  an  hour  of  averaging.  (Admittedly,  we  have  ignored  stratifica¬ 
tion  effects  on  these  estimates  for  simplicity  and  because  these  effects  are  not 
well  known  for  the  required  covariance  statistics.)  We  could,  of  course,  average 
the  fluxes  over  shorter  intervals  to  match  the  averaging  time  of  about  15  minutes 
we  estimated  for  the  wind  speed  and  temperature  gradients.  But  the  error  in  the 
fluxes  would  increase  with  shorter  averaging.  From  (3.5)  we  can  estimate  by  how 
much.  For  example,  reducing  the  averaging  time  fi'om  one  hour  to  15  minutes 
would  increase  the  flux  error  by  a  factor  of  2 — ^firom  10%  to  20%  in  our  analysis. 

We  have  undertaken  this  discussion  to  establish  the  averaging  constraints 
under  which  typical  Monin-Obukhov  similarity  fimctions  must  be  evaluated. 
Statistical  theory  and  experiment  suggest  that  it  may  be  possible  to  evaluate 
these  functions  from  only  15  minutes  of  averaging,  but  the  results  will  be  quite 
scattered  (cf  Haugen  et  al.  1971,  Wyngaard  1973,  p.  141),  especially  since  ^ 
includes  both  u,  and  t.  and  their  uncertainties  add  [see  (A3)].  In  fact,  we  know 
of  no  attempts  to  validate  Monin-Obukhov  similarity  theory  with  such  short 
averages.  Almost  all  published  similarity  fimctions  are  based  on  averages  of 
30-60  minutes. 

The  hope  is  that  scintillometer  measurements,  because  of  their  path-aver¬ 
aging,  might  provide  relief  from  these  averaging  constraints.  But  the  reality  is 
that  computing  surface  fluxes  from  scintillation  measurements  requires  using 
similarity  functions  derived  from  long  averaging  times.  No  evidence  exists  that 
Monin-Obukhov  similarity  theory  is  also  valid  for  the  1-10  minute  fluxes  often 
computed  from  scintillation  data  (e.g.,  Thiermatm  and  Grassl  1992,  De  Bruin  et 
al.  2002,  Hartogensis  et  al.  2002).  To  us,  assuming  that  Monin-Obukhov  simi¬ 
larity  functions  derived  from  30-60  minute  averages  are  equally  valid  when 
applied  to  short  averages  is  unjustified. 

If,  however,  the  path-averaging  in  scintillation  measurements  provides  such  a 
large,  fast  sample  that,  s^,  10  minutes  of  averaging  provides  flie  same  informa¬ 
tion  that  30-60  minutes  of  averaging  would,  using  existing  Monin-Obukhov 
similarity  functions  for  computing  scintillometer-derived  fluxes  may  be  justified. 
In  this  p^er,  we  thus  study  whether  short-term  samples  of  scintillometer  and 
£q  values  faithfully  represent  the  behavior  of  these  quantities  over  an  hour.  If  the 
answer  is  “yes”  in  general,  using  short  scintillometer  averages  and  the  Monin- 
Obukhov  similarity  fimctions  to  compute  short-term  fluxes  is  defensible  because 
the  path-averaging  has  mitigated  the  time-averaging  constraints  that  hamper  point 
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measurements.  If  the  answer  is  “no,”  however,  scintillometers  confer  no  time¬ 
averaging  benefits  because  their  data  must  still  pass  through  the  traditional 
Monin-Obukhov  similarity  fimctions  to  3deld  surface  fluxes. 
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4  AND  ^0  DISTRIBUTIONS 

To  fairly  compare  the  1997  and  1998  data  sets,  we  extr^olated  both  the  C„ 
and  ^5  values  to  a  reference  height  of  3  m  on  the  basis  of  (1.5),  (1.6),  and  (1.7). 
That  is. 


and 


(4.1) 


4(3)  =  4(z) 


^3.00 

<  z  J 


sl/d 


(4.2) 


where  (z)  and  f  q  (z)  are  the  values  measured  at  path  height  z,  and  (3) 
and  (3)  are  the  extrapolated  values  at  3  m.  Clearly,  (4.1)  and  (4.2)  ignore 
the  stability  dependencies  contained  in  g(C)  <|>e  (C)  •  But  since  the  extrap¬ 
olations  are,  at  most,  just  40  cm,  we  are  safe  with  the  approximations  that 
g(3/L)/g(z/L)  and  [(jis  (3/L)/(j)s  (z/L)]'^'*  are  bodi  near  1. 

Of  the  600  hours  of  scintillometer  data  that  survived  our  quality  controls, 
almost  all  of  the  data  collected  in  imstable  stratification  had  z/L  values  between 
-1  and  0,  with  the  majority  between  -0.1  and  0.  The  data  collected  in  stable 
stratification  featured  z/L  values  as  large  as  10,  but  most  z/L  values  were 
between  0  and  2. 

On  using  measurements  of  u.  and  Hs  from  another  Arctic  experiment, 
AIDJEX  (the  Arctic  Ice  Dynamics  Joint  Experiment),  Andreas  (1989a)  inferred 

from  a  relation  similar  to  (1 .6)  and  found  that  a  beta  distribution  fitted 
seasonal  histograms  of  these  values  well.  Figures  3-6  show  similar  seasonal 
histogram  of  our  SHEBA  €„  and  f  q  measurements  for  autumn  (namely, 
October  and  November  1997),  winter  (namely,  December  1997),  spring 
(namely.  May  1998),  and  summer  (June,  July,  and  August  1998). 

Briefly,  the  beta  probability  distribution  function  of  random  variable  y  is 
(e.g.,  Harr  1987,  p.  79  ff.) 


p(y) 


(b-a)  P(m  +  l,n  +  l) 


y-a 

b-a 


Nil. 


b-y 

b-a 


(4.3) 
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Here  a  and  b  are,  respectively,  the  lower  and  upper  bounds  on  y;  P  is  the  beta 
function  (Abramowitz  and  Stegun  1965,  p.  258), 


P(m  +  l,n  +  l) 


r(m+i)r(n+i) 
r(m  +  n  +  2) 


(4.4) 


and  r  is,  again,  the  gamma  function  (Abramowitz  and  Stegim  1965,  p.  255). 

On  calculating  the  sample  mean  (y )  and  sample  standard  deviation  (Sy)  of  y, 
we  define  two  new  variables 

y  =  (y-a)/(b-a)  (4.5) 

and 

V  =  [sy/(b-a)J  .  (4.6) 

These  then  predict  the  parameters  m  and  n: 


m  =  (y  /v)(l-y)-(l+y). 


(4.7) 


n  =  (m+l)/y  -  (m  +  2)  .  (4.8) 

In  Figures  3-6,  since  we  plotted  on  a  logarithmic  scale,  we  also  based 
our  fitted  histograms  on  the  logarithm  of  .  That  is,  in  calculating  the  sample 
mean  and  sample  standard  deviation  and  in  setting  the  a  and  b  values  in  die  beta 
distributions  for  the  seasonal  histograms,  we  used  hi(C^)  rather  than  C^. 

Because  of  instrument  problems  and  the  resulting  sparse  data  returns  in  some 
seasons,  our  seasonal  €„  histograms  in  Figures  3-6  aie  not  as  representative  of 
the  true  Arctic  seasons  as  Andreas’s  (1989a).  Rather,  we  offer  the  histograms 
in  Figures  3-6  as  further  evidence  that  the  beta  distribution  is  useful  for  repre¬ 
senting  die  variability  in  this  statistic.  The  £q  histograms  in  Figures  3-6  are,  we 
believe,  the  first  attempts  at  constructing  seasonal  histograms  far  £q. 
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Figure  3.  Histograms  of  the  minute-averaged  C^  and  values  measured 
at  SHEBA  during  autumn  (namely,  October  and  November  1997).  The 
number  in  the  upper  right  comer  gives  the  number  of  minute-averaged 
values  used  to  create  the  histogram.  The  solid  curve  is  a  beta  distribution 
fitted  to  the  histogram.  In  the  C^  plot,  the  lower  and  upper  limits  of  the 
beta  distribution,  a  and  b,  are  and 

respectively;  in  the  £g  plot,  these  limits  are  0  and  20  mm. 
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Figure  4.  As  in  Figure  3  except  for  winter  (that  is,  December  1997). 
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Figure  5.  As  in  Figure  3  except  for  spring  (that  is,  May  1998). 
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Figure  6.  As  in  Figure  3  except  for  summer  (June,  July,  August  1998). 
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Although  the  beta  distribution  requires  four  parameters  (i.e.,  y ,  Sy,  a,  b), 
rather  than  the  two  needed  for  normal  or  lognormal  distributions,  for  example, 
the  fit  is  not  especially  sensitive  to  a  and  b.  In  the  four  histograms  in  Figures 
3-6,  we  used  the  same  a  and  b  values.  In  (lO"'* )  and  In  (l0"‘^ ) ,  respectively. 
These  are  also  the  same  values  that  Andreas  (1989a)  used.  Likewise,  for  the  fom 
histograms,  we  used  Are  same  limits,  0  and  20  mm.  We  conclude  that  the  beta 
distribution  is  a  useful  tool  for  representing  broad  or  narrow  histograms  that 
result  from  few  or  many  data  samples.  And  for  our  Arctic  scintillometer  data, 
the  fit  is  fairly  insensitive  to  the  choices  of  a  and  b. 

Although  we  have  no  theoretical  justification  for  choosing  the  beta  distribu¬ 
tion,  Harr  (1987,  p.  77  ff.),  for  example,  provides  ample  practical  justification. 
Briefly,  with  different  choices  of  m  and  n,  the  beta  distribution  can  be  symmetric 
or  skewed;  it  can  approximate  a  normal  or  a  lognormal  distribution;  and  it  can 
even  represent  a  bathtub  distribution. 

The  lognormal  distribution  is  also  commonly  used  to  model  propagation 
statistics  (e.g.,  Ben-Yosef  and  Goldner  1988,  Hill  and  Frehlich  1997).  In  parti¬ 
cular,  Frehlich  (1992)  found  “satisfactory”  agreement  when  he  used  the  log¬ 
normal  distribution  to  represent  bodi  Cn  and  measurements  comprising  512 
consecutive  4-s  averages.  For  comparison,  we  thus  plot  autumn  and  spring 
histograms  again  in  Figures  7  and  8  but  now  fitted  with  lognormal  distributions. 

Since  the  axes  in  Figures  7  and  8  are  logarithmic,  the  plotted  lognormal 
distributions  are  nearly  S5mimetric.  (Lognormal  distributions  plotted  as  contin¬ 
uous  functions  of  €„  would  be  perfectly  symmetric;  the  need  for  bin-averaging 
in  these  plots  distorts  the  S}mimetry  a  bit.)  The  actual  histograms,  however, 
are  decidedly  skewed.  As  a  result,  in  both  figures,  the  lognormal  distribution 
imderestimates  the  location  of  the  peak  in  the  histogram,  undeipredicts  the  occur¬ 
rence  of  the  larger  €„  values,  and  has  an  upper  tail  with  too  much  area.  Andreas 
(1989a)  identified  similar  faults  during  comparisons  of  beta  and  lognormal  fits  of 
histograms.  Thus,  aside  from  computational  ease,  the  lognormal  distribution 
has  nothing  to  recommend  it  over  the  beta  distribution  for  fitting  histograms. 

Seasonal  histograms  of  €„  and  f  q  » however,  are  not  really  our  focus  here; 
we,  thus,  present  Figures  3-6  mainly  to  segue  from  Andreas’s  (1989a)  work  to 
ours  and  as  background  on  using  beta  distributions  for  representing  scintillation 
statistics.  Here  we  necessarily  concentrate  instead  on  the  short-term  behavior  of 
and  f  Q  and  find  the  beta  distribution  equally  useful. 
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Figure  7.  histogram  from  Figure  3,  except  here  the  soiid  curve  is  a  log- 
normal  distribution  fitted  to  the  histogram. 
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Figure  8.  histogram  from  Figure  5,  except  here  the  solid  curve  is  a  log¬ 
normal  distribution  fitted  to  the  histogram. 
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Figures  9-1 1  show  histograms  and  beta  distributions  fitted  to  only  one  hour 
(i.e.,  60  one-minute  averages)  of  and  values  for  three  different  occasions. 
Later  in  the  paper,  we  will  base  statistical  analyses  of  these  three  periods  on  beta 
distributions.  Therefore,  we  present  Figures  9-1 1  to  emphasize  that  beta  distribu¬ 
tions  can  fit  both  broad  and  narrow  histograms  comprising  far  fewer  data  points 
than  in  the  seasonal  histograms  in  Andreas  (1989a)  or  in  Figmes  3-6. 

To  make  Figures  9-1 1,  we  imposed  the  same  lower  and  upper  limits  on  both 
the  Cn  and  distributions  that  we  did  in  Figures  3-6.  Yet  even  with  these 
constant  limits,  the  beta  distribution  is  able  to  represent  both  wide  (e.g.,  Figmes 
3, 5,  and  11)  and  narrow  (e.g..  Figures  4, 9,  and  10)  histograms.  This  ability  of 
the  beta  distribution  with  fixed  lower  and  upper  limits  to  capture  the  features  of 
the  histogram,  even  when  there  are  few  points,  is  important  for  om  later  analysis. 
Also,  by  fixing  a  and  b  for  our  and  analyses,  we  have  essentially  tiuned 
the  beta  distribution  into  a  two-parameter  model. 

Another  important  feature  of  the  beta  distribution  for  our  application  is  that 
it  is  bounded  by  a  and  b.  The  winter  f  q  histogram  (i.e..  Figure  4)  looks  bell 
shaped;  we  could  have  tried  fitting  it  with  a  normal  distribution.  But  given  the 
sample  mean  and  sample  standard  deviation  of  these  winter  f  ^  measiuements, 
die  normal  distribution  would  have  predicted  a  finite  probability  that  could 
be  negative — a  physical  impossibility.  By  using  a  beta  distribution  for  q  and 
setting  the  a  value  to  0,  we  prevent  this  imphysical  result. 
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Figure  9.  Histograms  of  60  minute-averaged  and  values  measured  at 
SHEBA  between  07  and  08  UT  on  28  November  1997.  The  solid  curve  is  the 
beta  distribution  fitted  to  the  histogram.  In  the  plot,  the  lower  and  upper 
limits  of  the  beta  distribution,  a  and  b,  are  ln(10"^®  m"*'®)  and  ln(10"''*  m"^), 
respectively;  in  the  £q  plot,  these  limits  are  0  and  20  mm. 
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Figure  10.  As  in  Figure  9,  except  these  data  are  for  09-10  UT  on  28 
November  1997. 
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Figure  11.  As  in  Figure  9,  except  these  data  are  for  07-08  UT  on  23  May 
1998. 
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5  AVERAGING  AND 

With  the  knowledge  that  beta  distributions  can  reliably  represent  histograms 
of  Cn  and  ,  we  can  set  confidence  intervals  for  measurements  of  and  . 
In  other  words,  we  can  investigate  the  role  of  averaging  time  in  determining  the 
confidence  we  have  in  various  averages  of  and  .  Remember,  the  funda¬ 
mental  question  we  are  researching  is  whether  short-term  averages  of  and 
can  yield  meaningfiil  estimates  of  u,  and  H5.  Our  view,  again,  is  that,  because  of 
the  need  to  use  empirical  similarity  functions  to  compute  these  fluxes,  the  answer 
is  “yes”  only  when  the  short-term  measurements  of  and  faithfully  repre¬ 
sent  the  behavior  of  the  and  time  series  over  30-60  minutes. 

Define  a  time-dependent  average  X(t)  and  standard  deviation  s(t)  as 


=  Tip' 


and 


'(t)  = 


N-1 


-NX^(t) 

.  i=l 


(5.1) 


(5.2) 


Here  Xj  is  the  time  series  of  minute  averages  of  either  or  .f  q  ,  t  =  (N  -  l)At, 
and  At  is  one  minute  for  our  data.  That  is,  X(t)  and  s(t)  are  the  sample  mean  and 
sample  standard  deviation  for  all  data  in  a  time  series  from  the  first  point  to  the 
Nth  point,  which  corresponds  to  time  t. 

Figure  9  shows  histograms  for  the  and  £q  values  collected  on  28 
November  1997  between  07  and  08  UT.  Figure  12  shows  the  hour-long  time 
series  of  these  minute-averaged  and  £q  values.  The  figure  also  includes 

the  averages  (t)  and  £(,  (t)  computed  wifli  (5.1)  between  time  zero  and 
time  t  and  90%  confidence  intervals  on  the  and  £0  populations  at  each 
minute  based  on  beta  distributions  computed  fi'om  the  average  pC(t)]  and  the 
standard  deviation  [s(t)]  of  each  variable  at  time  t.  [As  before,  for  ,  all 
statistics  are  based  on  In  ) .]  The  lower  and  upper  limits  of  the  beta  dis¬ 
tribution  for  here  again  correspond  to  10”’®  and  10”’^  m"^®;  for 
limits  are  again  0  and  20  mm. 


Flux  Averaging 


27 


Time,  t  (minutes) 

Figure  12.  Hour-long  time  series  of  the  minute-averaged  values  of  €„  and 
^0  and  averages  of  these  computed  from  (5.1)  from  time  zero  until  the 
plotted  time,  t.  The  errors  bars  show  90%  confidence  intervals  based  on  a 
beta  distribution  computed  from  the  average  and  standard  deviation  at  time 
t  [i.e.,  from  (5.1)  and  (5.2)]. 
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In  Figure  12,  we  show  90%  confidence  intervals  for  the  and  popula¬ 
tions  (as  opposed  to  confidence  intervals  on  the  respective  means)  since  these 
intervals  represent  the  variability  in  the  path-averaged  atmosphere  as  observed 
by  the  scintillometer.  Remember,  our  objective  is  to  see  how  effective  path¬ 
averaging  instruments  are  in  smoothing  out  nonstationarity  or  intermittency  in 
the  atmosphere. 

Both  the  Cn  and  £q  time  series  in  Figure  12  support  the  contention  that 
statistically  meaningful  values  of  these  quantities  can  be  obtained  quite  quickly. 
Here,  after  only  10  minutes  of  averaging,  the  90%  confidence  intervals  for  both 
the  Cn  and  £q  populations  bracket  the  respective  averages  for  the  remainder  of 
the  hour  and  also  bracket  almost  all  of  the  subsequent  minute-averaged  values. 
That  is,  the  first  10  minutes  of  data  faithfully  represent  all  the  and  £q  values 
observed  diuing  the  hour.  For  such  series,  using  short-term  averages  of  and 
£q — say  10-minute  averages — ^to  make  inertial-dissipation  estimates  of  u»  and 
Hs  seems  reasonable.  Averaging  and  £q  for  30-60  minutes  to  match  the 
averaging  represented  in  the  Monin-Obukhov  similarity  functions  necessary 
for  the  inertial-dissipation  estimates  would  not  have  changed  the  C\  and  £q 
averages  significantly. 

To  confirm  this  conclusion,  we  have  estimated  u.  and  Hs  from  these  time 
series  using  (1.5)-(1.7)  and  the  similarity  functions  g(0  and  (j)^  (^)  given  in 
Appendix  A.  The  averages  of  the  first  10  minutes  of  and  f  q  values  give  u» 
and  Hs  values  of  0.0745  m  s~’  and  -6.68  W  m“^,  respectively,  while  the 
and  £^  averages  for  the  entire  60-mmute  series  yield  values  of  0.0754  m  s”’  and 
-6.78  W  m"^.  That  is,  the  differences  here  between  fluxes  based  on  60-minute 
averages  and  averages  of  only  the  first  10  minutes  of  the  series  are  about  1%  for 
both  u»  and  Hs.  In  fact,  the  u,  and  Hs  values  computed  from  the  and  £^ 
averages  for  each  10-minute  block  within  these  60-minute  series  are  within  about 
3%  of  the  u,  and  Hs  values  computed  from  the  hourly  averages  of  and  f  q  . 
Again,  these  results  support  the  premise  that  path-averaging  instruments  provide 
accmate  estimates  of  the  surface  fluxes  quickly. 

The  histograms  in  Figure  10  and  the  time  series  in  Figure  13,  however, 
present  a  different  situation  just  two  horns  later.  During  this  hour,  has  a 
downward  trend,  while  £q  has  an  upward  trend.  Both  series  also  show  shorter 
periods  of  nonstationarity.  As  a  result,  the  90%  confidence  interval  for  the 
population  calculated  from  10  minutes  of  data  brackets  fewer  than  half  of  the 
subsequent  minute  averages  and  does  not  even  include  the  average  at  the  end 
of  the  hour.  Likewise,  the  90%  confidence  interval  for  the  population  ob¬ 
tained  from  10  minutes  of  averaging  brackets  fewer  than  half  of  the  subsequent 
minute  averages  of  f  (,  and  barely  includes  the  f  q  average  at  the  end  of  the  hour. 
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In  other  words,  for  the  data  in  Figures  10  and  13, 10  minutes  of  averaging  does 
notfaiihfully  predict  the  variability  of  the  and  values  that  occur  during 
this  hour. 
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Figure  13.  As  in  Figure  12  but  two  hours  later. 
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This  example  thus  refutes  the  hypothesis  that  short-term  averages  of  Q\  and 
f  0  can,  in  general,  yield  meaningful  estimates  of  the  surface  momentiun  and  sen¬ 
sible  heat  fluxes.  Here,  because  of  the  nonstationarity,  it  is  not  obvious  that  short¬ 
term  Cn  and  f  0  averages  can  be  meaningfully  convolved  with  functions  based 
on  30-60  minutes  of  averaging.  In  essence,  nonstationarity  can  confound  path¬ 
averaging  sensors  just  as  it  confounds  point  sensors. 

To  demonstrate,  we  again  compare  u,  and  Hs  values  computed  from  the 
first  10  minutes  of  these  and  series  and  from  the  entire  series.  and 
f  Q  averages  for  the  first  10  minutes  yield  u*  and  H  values  of  0.0844  m  s"’ 
and  -8.18  W  m"^,  respectively,  while  the  60-minute  averages  yield  values  of 
0.0795  m  s”'  and  -7.18  W  m"^.  These  u»  values  differ  by  6.1%,  while  the  Hs 
values  differ  by  13.9%.  The  fourth  10-minute  block  in  these  series  produced 
worse  results;  the  u,  and  Hs  values  based  on  these  10-minute  averages  differ 
from  values  based  on  the  60-minute  averages  by  7.8%  and  15.6%,  respectively. 

While  these  differences  in  the  Hs  estimate,  in  particular,  are  not  large  in  mag¬ 
nitude,  the  percentage  differences  are  significant.  In  mid-latitudes,  over  land, 
where  Hs  can  be  200-300  W  m'^,  for  example,  a  15%  rmcertainty  in  Hs  amoxmts 
to  30-45  W  m"^.  Such  an  uncertainty  is  much  too  large  if  we  have  intentions,  for 
instance,  to  ever  use  scintillation  to  improve  climatological  estimates  of  the  sur¬ 
face  heat  budget,  where  an  accuracy  of  5-10  W  m"^  is  essential  (e.g.,  Mitchell 
1989,  Kiehl  and  Trenberth  1997). 

Figure  14  shows  another  horn  of  €„  and  values  from  May  1998  to 
emphasize  these  points.  (The  corresponding  histograms  are  in  Figure  11.)  In 
contrast  to  Figure  13,  the  and  time  series  here  do  not  have  trends  but  do 
display  a  lot  of  variability.  As  a  result,  although  the  confidence  interval  on  the 
cl  population  at  10  minutes  brackets  all  subsequent  averages,  it  brackets 
only  about  half  of  the  subsequent  minute  averages.  The  situation  is  not  quite  as 
bad  for  the  £q  series,  mainly  because  variability  early  in  the  series  leads  to  a 
large  90%  confidence  interval  for  the  population  after  10  minutes  of  averaging. 
Still,  20%  of  the  subsequent  minute-averaged  values  are  outside  this  90% 

confidence  interval.  And  the  confidence  interval  for  the  population  after  an  hour 
of  averaging  is  25%  larger  than  the  confidence  interval  at  10  minutes. 

Hence,  again  10  minutes  of  scintillometer  data  did  not  provide  a  very 
accurate  picture  of  the  and  £q  populations  we  encountered  during  fliis  hour. 
Therefore,  it  does  not  seem  prudent  to  use,  say,  the  10-minute-averaged  €„  and 
f  Q  values  from  these  series  in  inertial-dissipation  estimates  of  the  smface  fluxes. 
The  nonstationarity  leads  to  confidence  intervals  for  the  and  £q  populations 
with  widths  that  are  significant  fractions  of  the  respective  mean  values:  about 
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To  confinn  this  conclusion,  we  again  compare  calculations  of  u,  and  Hs 
based  on  the  first  10  minutes  of  these  €„  and  times  series  and  on  the  entire 
60-minute  series.  The  first  10  minutes  yields  u«  and  H  values  of  0.0451  m  s“* 
and -1.32  Wm"^  respectively,  while  the  60-minute  and  averages  yield 
values  of  0.0488  m  s”'  and  -1.42  W  m  These  values  differ  by  7-8%.  Averages 
of  and  Iq  from  the  second  10-minute  block  in  these  series,  however,  yield 
u»  and  Hs  values  that  differ  from  the  60-minute  averages  by  22%  and  24%, 
respectively.  These  are  large  exclusions  from  the  long-term  means  despite  the 
fact  that  these  data  from  10-19  minutes  were  also  used  to  compute  the  60-minute 
averages. 

The  idea  that  path-averaging  instruments  can  yield  meaningful  estimates 
of  turbulent  surface  fluxes  in  a  fraction  of  the  time  that  point  measurements  can 
springs  from  the  hope  that  the  path-averaging  would  quickly  sample  enough 
turbulent  eddies  to  minimize  the  effects  of  nonstationarity.  As  Figures  13  and  14 
demonstrate,  for  propagation  paths  of  300-350  m  at  least,  this  does  not  happen  in 
general.  Nonstationarity  is  still  a  problem  for  path-averaging  instruments. 
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6  QUANTIFYING  THE  NONSTATIONARITY 

Since  nonstationarity  turns  out  to  present  sampling  problems  for  path¬ 
averaging  sensors  as  well  as  for  point  sensors,  we  use  our  scintillometer  data 
to  study  w^s  to  quantify  this  nonstationarity.  Mahrt  (1998)  defines  a  “non¬ 
stationarity  ratio”  just  for  this  purpose.  We  adapt  his  method  to  our  data. 

Again,  let  x  represent  the  minute-averaged  time  series  of  or  values. 
Consider  an  horn  series  that,  thus,  includes  60  of  these  values.  Block  these  60 
values  into  I  records  of  J  samples  each.  The  x  series  thus  now  has  two  subscripts, 
i  and  j  (i.e.,  Xy),  where  i  goes  from  1  to  I  and  j  goes  from  1  to  J. 

Mahrt  (1998)  defmes  the  ‘Svithin-record”  standard  deviation  of  x  for  record  i 
as 


where 


Xi 


(6.1) 


(6.2) 


He  then  averages  Swi.i  over  the  I  records  to  obtain  one  estimate  of  the  within- 
record  standard  deviation. 


(6.3) 


From  this  statistic,  he  defmes  the  random  error  as 


RE  =  s„./J’ 


1/2 


(6.4) 


Mahrt  (1998)  next  defmes  the  “between-record”  standard  deviation  as 


(6.5) 
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where 


1  i=l  IJ  i=l  j=l 


(6.6) 


Clearly,  in  our  analysis,  I J  =  60.  Mahrt’s  nonstationarity  ratio  is  then  just 


■t\n 

nr  =  ^  =  (6.7) 

RE  s 

*  ‘^Wl 

Conceptually,  Sbtw  is  large  when  nonstationarity  produces  long  segments  of 
the  series  that  remain  above  or  below  the  series  mean,  as  in  both  the  and 
series  in  Figure  13.  On  the  other  hand,  Swi  is  large  and  NR  tends  to  be  smaller 
when  the  time  series  has  a  great  deal  of  random  variability,  as  in  the  €„  plot  in 
Figme  14.  Thus,  Mahrt’s  (1998)  nonstationarity  ratio  compares  the  effects  of 
coherent  behavior  with  the  inherent  variability  in  the  time  series.  According  to 
Mahrt,  “for  stationary  conditions,  NR  is  approximately  unity.”  NR  values  signi¬ 
ficantly  larger  than  one  reflect  nonstationarity. 

We  have  calculated  NR  twice  each  for  the  raw  60-minute  and  time 
series  plotted  in  Figures  12-14  to  see  if  this  ratio  can  complement  our  analysis  in 
the  previous  section.  In  one  case,  1  =  6  and  J  =  10;  and  in  the  second  case,  I  =  10 
and  J  =  6.  Table  1  summarizes  the  results. 

The  first  thing  we  notice  in  the  table  is  fliat  the  NR  values  calculated  for  the 
two  I-J  pairings  are  not  necessarily  the  same.  For  the  and  £q  traces  in 

Figures  12  and  14,  the  respective  NRs  for  the  two  I-J  pairs  are  similar,  if  not 
identical.  But  for  Figure  13,  the  NRs  for  ,  7.05  and  5.27,  and  for  f  q,  6.54  and 
5.36,  each  differ  by  over  20%.  Mahrt  (1998)  does  not  suggest  any  obvious  guide¬ 
lines  for  choosing  I  or  J;  but  from  the  basis  for  his  method,  we  believe  that  J 
should  represent  a  period  that  is  typical  of  the  large-scale  variability  in  the  signal. 
The  and  traces  in  Figure  13  have  large-scale  coherence  for  periods  longer 
than  6  minutes.  As  a  result,  the  J  =  10  and  J  =  6  calculations  yield  significantly 
different  results  for  these  series. 

On  the  other  hand,  the  nonstationarity  ratios  listed  in  Table  1  both  confirm 
some  of  our  previous  conclusions  that  we  based  on  confidence  intervals  and 
contrast  with  them.  For  example,  the  and  £q  traces  in  Figure  12  are  well 
behaved,  and  we  had  concluded  that  short-term  averages  from  these  series  would 
be  useful  for  estimating  fluxes.  Likewise,  tire  NR  values  for  these  series  for  both 
I-J  pairs  are  consistent  and  are  two  of  the  three  smallest  collections  in  tire  table. 
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In  particular,  the  NRs  for  the  £q  series,  1.6-1. 7,  are  near  one,  the  value  Mahrt 
(1998)  associates  with  a  stationary  time  series. 


Table  1.  Calculations  of  the  nonstationarity  ratio,  NR,  from  (6.7)  for  the 
raw  and  series  shown  in  Figures  12-14  for  two  choices  of  I  and  J. 
“Average”  denotes  the  one-hour  average.  “Crossings”  counts  the  num¬ 
ber  of  times  the  series  crosses  this  average.  “59/Crossings”  is  another 
measure  of  nonstationarity  explained  in  the  text. 

28  Nov  1997 

28  Nov  1997 

23  May  1998 

07-08  UT 

(Figure  12) 

09-10  UT 

(Figure  13) 

07-08  UT 

(Figure  14) 

Average 

4.14x10-’“ 

4.01x10-’“ 

3.27x10-’* 

NR 

1  =  6,  J  =  10 

2.75 

7.05 

2.23 

1  =  10,  J  =  6 

2.97 

5.27 

2.23 

Crossings 

20 

9 

18 

59/Crossings 

2.95 

6.55 

3.28 

^0 

Average  (mm) 

6.62 

6.45 

10.54 

NR 

I=6.  J  =  10 

1,63 

6.54 

3.77 

1  =  10,  J  =6 

1.68 

5.36 

3.35 

Crossings 

22 

7 

13 

59/Crossings 

2.68 

8.43 

4.54 

The  nonstationarity  ratios  also  corroborate  our  evaluation  of  the  and  £q 
series  in  Figure  13.  These  both  had  unacceptable  short-term  behavior,  largely 
because  of  their  trends;  and  the  NRs  for  these  are  the  four  largest  values  in  the 
table.  Evidently,  if  a  series  is  stationary,  the  choice  of  I  and  J  is  not  too  crucial; 
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Both  I-J  pairs  yielded  comparable  NR  values  for  both  the  and  traces  in 
Figure  12.  But  when  a  series  is  nonstationary,  as  is  Figure  13,  different  choices 
for  I  and  J  can  yield  markedly  different  values  of  NR  These  NR  values, 
nevertheless,  are  large  and  do  suggest  what  they  should — ^that  the  series  is 
nonstationary. 

Lastly,  the  nonstationarity  ratios  for  the  and  series  in  Figure  14 
suggest  conclusions  that  contrast  with  our  earlier  discussion  of  these  series.  We 
had  also  judged  these  two  series  uimeliable  for  producing  short-term  averages 
because  of  their  large  variability.  The  NRs  of  the  series  in  Figure  14, 
however,  both  2.23,  are  smaller  than  the  NR  of  the  supposedly  “good”  series 
in  Figure  12.  And  the  NR  values  for  the  series,  3.77  and  3.35,  are  larger  than 
the  NRs  for  the  €„  series,  although  we  had  earlier  concluded  fliat  the  series 
in  Figure  14  was  “better”  for  averaging  than  the  series  on  the  basis  of  oiu: 
uncertainty  analysis.  These  discrepancies  raise  questions  about  what  nonstation¬ 
arity  means  in  general  and  what  Mahrt’s  (1998)  nonstationarity  ratio  quantifies  in 
particular. 

In  reviewing  (6.1)-(6.7),  we  see  that  NR,  through  its  linear  dependence  on 
the  between-record  variability,  Sbtw,  tends  to  be  large  if  coherent  sections  of  the 
record  are  above  or  below  the  mean.  On  the  other  hand,  because  of  its  inverse 
dependence  on  within-record  variability,  Swi,  NR  tends  to  be  small  if  the  time 
series  has  large-amplitude  random  variability.  The  trace  in  Figure  14  has 
such  large-amplitude  variability  (i.e.,  large  Swi),  but  the  events  are  not  necessarily 
coherent  (i.e.,  small  Sbtw).  As  a  result,  the  NR  for  this  series  is  relatively  small 
despite  error  bars  that  increase  wifli  time.  The  trace  in  Figure  14,  in  contrast, 
has  large-amplitude  variability  (i.e.,  large  Swi),  but  that  variability  is  more  coher¬ 
ent.  There  are,  essentially,  four  segments  in  the  time  series:  two  above  and  two 
below  the  mean.  Thus,  Sbtw  is  also  large,  and  this  combination  yields  mid-range 
NRs  that  are  larger  than  the  NRs  for  the  series. 

In  summary,  according  to  Mahrt’s  (1998)  definition  of  nonstationarity,  as 
formalized  in  his  nonstationarity  ratio  NR,  two  effects  contribute  most  to  non¬ 
stationarity:  a  trend,  and  large,  coherent  excursions  fi'om  the  mean.  By  “large” 
here,  we  mean  events  that  span  a  significant  fraction  of  the  time  series,  say  10% 
of  it.  We  are  tempted  to  call  such  events  “intermittency,”  but  the  definition  of  this 
term  in  the  new  Glossary  of  Meteorology  (Glickman  2000)  does  not  seem  to  fit 
our  description.  Consequently,  we  will  simply  call  these  nonstationary  events. 

Again  according  to  Mahrt  (1998),  large-amplitude,  random  variability  is  the 
opposite  of  nonstationarity:  This  variability  actually  reduces  his  nonstationarity 
ratio.  In  essence,  if  a  series  is  highly  variable,  coherent  exclusions  from  the  mean 
should  not  be  surprising  and,  thus,  should  not  be  judged  as  nonstationarity.  In 
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our  luicertainty-based  analysis,  however,  such  highly  variable  series  lead  to  wide 
confidence  intervals  and,  in  turn,  to  uncertain  estimates  of  the  surface  fluxes.  In 
other  words,  nonstationarily  (at  least  as  defined  by  Mahrt)  is  not  the  sole  cause  of 
uncertain  flux  estimates:  Random  variability  also  leads  to  uncertainty. 

These  insights  into  Mahrt’s  (1998)  nonstationarity  ratio  led  us  to  realize  that, 
generally,  a  stationary  time  series  exhibits  many  zero-crossings,  while  a  non¬ 
stationary  series  has  long  segments  above  or  below  zero.  (Here,  zero-crossings 
refer  to  a  series  with  the  mean  removed.)  Table  1  also  lists  the  number  of  zero- 
crossings  for  the  raw  €„  and  £q  series  in  Figures  12-14,  where  that  count  in¬ 
cludes  both  up-going  and  down-going  zero-crossings.  These  coxmts  mirror  our 
analysis  based  on  NR.  The  and  £q  series  in  Figure  12  have  many  crossings 
and  small  NR;  the  series  in  Figure  13  have  few  crossings  and  large  NR;  and  the 
series  in  Figure  14  have  intermediate  numbers  of  crossings  and,  generally,  inter¬ 
mediate  NR. 

A  series  of  60  samples  can  cross  zero  at  most  59  times  and  must  cross  zero 
at  least  once.  Denote  the  number  of  crossings  as  C  and  this  maximum  number 
of  crossings  as  M.  Clearly,  1<M/C<M .  Mahrt’s  (1998)  NR,  in  contrast,  is  not 
necessarily  always  greater  than  or  equal  to  one.  Although  Mahrt  states  fliat  his 
NR  is  “approximately”  one  for  a  stationary  series,  we  have  created  artificial 
series  with  NR  less  than  one  (for  example,  a  60-point  square  wave  of  amplitude 
1  and  wavelength  12).  It  seems  useful  to  have  a  nonstationarity  metric  with 
obvious  lower  and  upper  limits.  In  counting  zero-crossings,  M/C  =  1  would 
mean  that  each  consecutive  point  switches  from  above  to  below  the  mean  or  vice 
versa.  This  behavior  is  a  good  indication  that  the  series  is  stationary.  In  contrast, 
M/C  =  M  would  mean  the  series  crosses  zero  only  once  and,  thus,  has  a  trend. 

The  M/C  values  in  Table  1  (i.e.,  “59/Crossings”)  are  similar  in  both  mag¬ 
nitude  and  trend  to  the  corresponding  NRs  and,  therefore,  seem  to  be  a  non¬ 
stationarity  metric  that  is  as  good  as  NR.  Counting  zero-crossings  is  also  much 
easier  than  making  the  calculations  necessary  to  evaluate  NR.  Another  benefit 
is  that  counting  zero-crossings  obviates  the  need  to  make  arbitrary  choices  of  I 
and  J. 

To  develop  a  better  sense  of  what  this  zero-crossing  metric  tells  us  about  a 
time  series,  we  computed  M/C  for  hourly  intervals  of  the  €„  and  £q  time  series 
in  Figme  2.  Figme  15  shows  the  results. 


38 


ERDC/CRRELTR-03-24 


Figure  IS.  The  nonstationarity  metric  M/C  for  hourly  segments  of  the 
and  time  series  in  Figure  2,  where  M  is  always  59  (=  60  - 1)  and  C  counts 
the  zero-crossings  during  the  hour.  Values  are  plotted  at  the  mid-point  of 
the  hour.  Because  during  the  last  two  hours  the  series  in  Figure  2  had 
gaps,  we  calculated  no  M/C  values  for  these  hours. 

The  first  thing  this  figure  suggests  is  that,  on  this  d^,  was  much  more 
nonstationary  than  £q  was.  In  fact,  between  02  and  04  UT,  when  climbed 
steadily  dirough  two  orders  of  magnitude  (see  Figure  2),  its  M/C  value  was  at  its 
maximum,  59.  That  is,  had  only  one  zero-crossing  each  hour.  During  the 
steady  decline  from  16  to  19  UT  in  Figure  2,  also  had  large  M/C  values, 
almost  20  (i.e.,  3  zero-crossings  during  each  horn). 

The  laigest  M/C  value  for  the  f  ^  record  in  Figure  2  also  occurs  for  02-03 
UT,  when  £q  is  climbing  out  of  its  local  minimum.  The  £q  record  also  has  large 
M/C  values  for  00-01  UT  and  10-1 1  UT,  when  £q  undergoes  some  large,  coher¬ 
ent  excursions.  We  conclude  from  Figure  15  that  M/C  is  a  good  indicator  of  a 
trend  in  a  time  series  and  of  coherent  excursions  on  either  side  of  the  mean. 

Finally,  Figure  15  suggests  that  and  £q  were  both  stationary  during 
the  five-hour  period  1 1-16  UT.  Here,  the  M/C  values  for  both  and  £q  are 
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between  2  and  6.  Likewise,  the  and  £(,  traces  in  Figure  2  show  fairly  con¬ 
stant  levels  here  and  more  modest  fluctuations  than  in  the  rest  of  the  records.  We 
thus  recommend  considering  a  time  series  to  be  stationary  if  M/C  is  less  than  6. 
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7  CONCLUSIONS 

Our  SHEBA  scintillometer  data  confirm  Andreas’s  (1989a)  conclusion 
that  a  beta  distribution  with  lower  and  upper  limits  corresponding  to  10"’®  and 
10"’^  m“^  reliably  models  the  distribution  of  Arctic  near-surface  values. 

We  are,  evidently,  the  first  to  report  that  a  beta  distribution  with  lower  and  upper 
limits  of  0  and  20  mm  is  also  useful  for  representing  values  of  the  inner  scale  of 
turbulence, 

Using  this  knowledge  of  the  sampling  distributions  for  €„  and  f  o » 
investigated  the  question  of  whether  short-term  averages  of  path-averaged 
and  f  0  values  can  yield  reliable  estimates  of  the  smface  stress  and  the  surface 
sensible  heat  flux  through  an  inertial-dissipation  calculation.  The  answer  turns  on 
the  fact  that  the  similarity  functions  used  in  the  inertial-dissipation  calculations 
derive  firom  point  measurements  based  on  averages  of 30-60  minutes.  The 
cautious  approach  is  to  convolve  short-term  averages  of  and  (i.e.,  on  the 
order  of  minutes)  with  these  similarity  functions  only  if  short  records  of  and 

f  Q  fairly  represent  the  behavior  of  time  series  over  30-60  minutes. 

Our  analyses  of  three  typical  hour-long  time  series  of  and  f  q  measme- 
ments  from  SHEBA  imply  that  path-averaging  confers  no  general  benefits  when 
the  objective  is  to  shorten  the  averaging  time  for  flux  sampling.  Short-term  series 
simply  do  not  alwtys  bracket  the  behavior  of  30-60  minute  series.  If  the  and 

f  0  time  series,  however,  are  stationary  over,  say,  hourly  intervals,  the  fluxes  of 
heat  and  momentum  computed  over  10-minute  subrecords  can  be  statistically 
reliable.  That  is,  for  stationary  series,  it  is  safe  to  convolve  short-term  and  f  q 
averages  with  similarity  functions  derived  fi'om  longer  averaging. 

But  nonstationarity  leads  to  uncertain  flux  measurements  with  path-aver¬ 
aging  instruments  just  as  it  does  for  point  measurements,  at  least  for  averaging 
paths  of 300-350  m.  On  the  basis  of  our  three  examples  and  the  averaging 
estimates  we  reviewed  in  Section  3,  we  expect  estimates  of  x  and  Hs  to  easily 
vary  by  20%  or  more  between  adjacent  10-minute  subrecords.  This  is  the  random 
error  resulting  from  variability  in  and  f  q  .  Estimates  of  x  and  Hs  could  also 
suffer  from  bias  errors  caused  by  using  similarity  functions  improven  for  short¬ 
term  averaging.  Admittedly  diough,  in  some  ^plications,  both  errors  may  be 
tolerable. 

In  our  view,  the  best  way  to  ensure  reliable  short-term  flux  estimates  from 
scintillometer  data  is  to  first  validate  the  Monin-Obukhov  similarity  functions 
fi'om  measurements  of  short-term  averages.  That  is,  make  point  measurements  of. 
Sty,  10-minute  averages  of  the  vertical  wind  speed  and  temperature  gradients  and 
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the  fluxes  uw  and  wt  and  see  how  well  these  data  fit  known  fiinctions  for 
(jig  (^)  and  g(Q.  Though  the  computed  functions  will  be  quite  scattered,  perhaps 
we  will  see  the  same  (|)g  (^)  and  g(Q  functions  emerging  that  longer  averaging 
has  yielded.  If  so,  the  experiment  establishes  the  validity  of  Monin-Obukhov 
similarity  for  short-term  averages  and  puts  scintillometer  flux  analysis  on  firm 
theoretical  and  experimental  footing.  To  our  knowledge,  however,  this  experi¬ 
ment  has  never  been  done,  and  there  are  no  guarantees  that  it  will  work.  As 
Wyngaard  (1973,  p.  41)  emphasizes  and  as  the  discussion  in  Section  3  demon¬ 
strates,  “Nature,  not  man,  determines  what  is  short-term.” 

Because  nonstationarity  turned  out  to  be  one  key  for  deciding  whether  a  flux 
estimate  was  reliable,  we  evaluated  two  metrics  for  quantifying  nonstationarity. 
Mahrt’s  (1998)  nonstationarity  ratio  tended  to  corroborate  our  analysis  of  the 
utility  of  and  for  estimating  fluxes  fliat  we  based  on  assigning  confidence 
intervals.  That  is,  acceptable  and  series  generally  had  small  values  of  NR, 
while  imacceptable  series  had  larger  NRs.  We  did  identify  some  examples  that 
went  against  intuition,  however.  These  series  had  large-amplitude,  random  vari¬ 
ability  and,  thus,  produced  averages  with  large  uncertainty;  still,  their  NRs  were 
fairly  small.  One  manifestation  of  nonstationarity  is,  thus,  extended  segments  of 
the  series  either  above  or  below  the  mean.  High  random  variability,  on  the  other 
hand,  is  the  opposite  of  nonstationarity  but  still  degrades  flux  estimates. 

In  light  of  these  insights,  we  realized  that  the  number  of  zero-crossings,  C,  in 
a  series  is  another  useful  measure  of  nonstationarity.  A  benefit  of  this  metric  is 
that  a  series  of  M+1  data  points  can  display  no  more  than  M  zero-crossings  and 
must  have  at  least  one  zero-crossing.  Gjnsequently,  M/C  can  have  values  only 
from  1  to  M.  When  it  is  1,  consecutive  points  in  the  series  just  switch  back  and 
forth  between  positive  and  negative — an  obviously  stationary  series.  When  it  is 
M,  the  series  crosses  zero  only  once  because  of  a  strong  trend.  Mahrt’s  (1998) 
nonstationarity  ratio,  on  the  other  hand,  does  not  have  obvious  limiting  values, 
depends  sensitively  on  choices  of  die  blocking  values  I  and  J,  and  requires  more 
elaborate  calculations  than  simply  counting  zero-crossings.  Finally,  in  our  six 
examples,  the  M/C  values  are  remarkably  near  the  corresponding  NRs.  M/C  is 
thus  another  simple  measure  of  nonstationarity. 
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APPENDIX  A:  THE  SIMILARITY  FUNCTIONS  g(g  AND  UQ 

Hill  (1997)  provides  an  authoritative  review  of  the  algorithms  and  similarity 
fimctions  necessary  to  compute  u*  and  t*  from  scintillation  statistics.  Although 
the  similarity  functions  g(0  and  (|)g  (^)  that  we  use  in  (1.6)  and  (1.7)  are  not 
necessarily  his  “preferences  ”  no  data  set  has  yet  proved  sufficient  to  distinguish 
significant  differences  between  his  choices  and  ours. 

For  (|)g,weuse 


(|),(C)  =  [l  +  0.46(-CrJ''  (Al) 

in  unstable  stratification  (  ^  <  0  )  and 

i(C)  =  [l  +  2.3C"=]’'’  (A2) 

in  stable  stratification  (Q  >0).  These  functions  are  from  Wyngaard  and  Cote 
(1971),  but  Andreas  (1988b)  modified  their  original  coefficients  to  reflect  a  von 
Karman  constant  of  0.40  rather  than  the  0.35  that  assiuned.  Here  also, 


«  =  f  -  ^-T  .  (A3) 

L  T,  u: 

where  L  is  the  Obukhov  length,  g  is  the  acceleration  of  gravity,  and  Ty  is  the 
average  virtual  temperature. 

Stable  stratification  is  a  special  concern  of  ours  in  this  SHEBA  data  set. 
Andreas  (1998, 2002)  and  Jordan  et  al.  (1999)  discuss  the  Monin-Obukhov 
similarity  frmctions  for  stable  stratification  and  recommend  the  so-called  “Dutch” 
formulation  (Holtslag  and  De  Bruin  1988)  for  die  dimensionless  temperature 
profile  function, 

(j.^  (0  =  1  +  0.7C  +  0.75C(6  -  0.35Oexp(-0.35C) .  (A4) 

This  has  the  best  properties  among  the  available  functions  in  very  stable 
stratification  (Andreas  2002).  We  can  incorporate  the  desirable  properties  of  this 
function  into  our  g(Q  equation  by  using  the  expression  for  g(Q  that  Andreas 
(1988b)  derives  from  the  scalar  variance  budget.  That  is. 
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5-92(1.,  (g 

<i>r(g  ’ 

where,  for  stable  stratification,  we  use  (A2)  for  (|)j,  and  (A4)  for  (j), . 
In  unstable  stratification,  we  use  (Edson  and  Fairall  1998) 

g(C)  =  5.92(1 -8cr^ 


(A5) 


(A6) 


because  it  matches  (A5)  best  at  C  =  0. 
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